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OPHTHALMIC DATA MEASURING APPARATUS, OPHTHALMIC DATA 
MEASUREMENT PROGRAM AND EYE CHARACTERISTIC MEASURING APPARATUS 

Technical Field 
5 The present invention relates to an ophthalmic data 

measuring apparatus for measuring appropriate correction data 
and/or estimating visual acuity in daily life, an ophthalmic 
data measurement progreun and an eye characteristic measuring 
apparatus . 

10 

Background Art 

Conventionally, as a technique for measuring ocular 
correction data, measurement of S (Sphere), C (Cylinder) and 
A (cucis) by a ref ractometer has been carried out. Besides, 

15 recently, an eye characteristic measuring apparatus capable of 
measuring higher order aberrations has also been developed, and 
not only S, C and A on a line like, for example, a ring of <|)3 
mm as in a ref ractometer , but also S, C and A on a plane when 
a pupil diameter is made various sizes can be calculated from 

20 lower order aberrations. By the eye characteristic measuring 
apparatus like this, especially after a refraction correcting 
surgical operation or in an eye disease, values closer to 
prescription values of eyeglasses or contact lenses than the 
ref ractometer can be calculated (for example, see 

25 JP-A-2002-204785, JP-A-2002-209854 , JP-A-2002-306416 , 

JP-A-2002-306417, etc.). 

Besides, as an apparatus for displaying the visibility 
of a subjective eye at the time of correction or by the eye, 
an apparatus by the present applicant Is disclosed (for example, 

30 see JP-A-2001-120504 , JP-A-7-100107 ) . In these apparatuses, 
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for example, the visibility of a predetermined index is 
displayed on display means based on the measured optical 
characteristic of the eye to be examined (eye to be measured) . 



5 Disclosure of the Invention 

However, in the objective calculation results of the 
conventional eye characteristic measuring apparatus and 
prescription values of eyeglasses, contacts, or the like, there 
is a case where a difference from an appropriate value occurs, 

10 and there has been a case where they are insufficient as 
evaluation of S, C and A. Besides, conventionally, since a 
measurement is made using a fixed value as the pupil diameter 
of an eye to be examined, there has been a case where an 
appropriate prescription value corresponding to the pupil 

15 diameter of the eye to be excunined can not be obtained • 

Besides, in the conventional measurement, although the 
visibility of an index or the like is displayed, the visual 
acuity of the subjective eye is not estimated. Further, the 
prediction of the visibility is often the visibility under a 

20 generally used visual acuity measurement condition, and the 
visibility and the visual acuity under an environment of the 
eye to be examined, for exsunple, in daily life, have not been 
obtained. 

Besides, conventionally, when PSF or MTF on the retina 
25 is simply evaluated, there has been a case where it is very 
difficult to obtain an appropriate evaluation, that is, a value 
close to a subjective test. 

In view of the above, the invention has an object to 
calculate an optical characteristic corresponding to a pupil 
30 diameter of an eye to be examined and correction data close to 
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an optimum prescription value and to perform more accurate 
measurement . 

Besides, according to an object of the invention, in the 
measurement results of an eye characteristic measuring 
5 apparatus which can measure higher order aberrations, in the 
case where a higher order aberration is included, a lower order 
aberration corresponding to the time of objective complete 
correction is not made compensation correction data, optical 
performance is evaluated with, for example, a Strehl ratio or 

10 a phase shift, a lower order aberration amount by which the 
Strehl ratio becomes large and/or the phase shift becomes small 
is calculated, and compensation correction data of S, C, A and 
the like at that time is obtained, so that correction data close 
to the optimum prescription value of eyeglasses /contacts is 

15 obtained. 

Further, an object is to obtain correction data close to 
a subjective value by performing simulation of visibility of 
an index for eye examination to obtain an appropriate correction 
element • 

20 The invention has an object to estimate visual acuity of 

an eye to be examined in luminance corresponding to an 
environment of a subjective eye in daily life (for example, in 
the daytime or in a room) . Besides , the invention has an object 
to estimate visual acuity with respect to an index of high 

25 contrast and/or low contrast in view of a pupil diameter of an 
eye to be examined in daily life. The invention has an object 
to predict contrast sensitivity in view of a pupil diameter. 
Besides, the invention has also an object according to which 
a pupil dicuneter in luminance corresponding to an environment 

30 of a subjective eye is used, correction data close to an optimum 



prescription value under the environment is obtained, and the 
visual acuity under the environment of the subjective eye at 
the time of correction by the obtained correction data is 
estimated. Besides, the simulation of an index, such as a 
5 Landolt's ring, on the retina in view of the size of a pupil 
area calculated in the middle of the process is also singly 
effective. 

According to the first solving means of this invention, 
there is provided, an ophthalmic data measuring apparatus 
10 comprising: 

a first light source part to emit a light flux of a first 
wavelength; 

a first illiuninatlng optical system for performing 
Illumination to condense the light flux from the first light 

15 source part on a vicinity of a retina of an eye to be examined; 

a first light receiving optical system for receiving a 
part of the light flux reflected by and returning from the retina 
of the eye to be examined through a first conversion member to 
convert it into at least substantially 17 beams; 

20 a first light receiving part for receiving the received 

light flux of the first light receiving optical system; and 
a calculation section to perform Zernike analysis based 
on an inclination angle of the light flux obtained by the first 
light receiving part, to obtain an optical characteristic of 

25 the eye to be examined, and (1) to estimate one of or two or 
more of a visual acuity, the optical characteristic and a 
sensitivity of the eye to be examined under an observation 
condition corresponding to an environment of the eye to be 
excunined, or (2) to calculate appropriate correction data 

30 suitable for the eye to be examined, 
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wherein the calculation section comprises: 
first means for obtaining measurement data indicating a 
refractive power distribution of the eye to be examined and 
pupil data including a value of a pupil diameter of the eye to 
5 be examined or a pupil diameter image and for obtaining lower 
order aberrations and higher order aberrations based on an 
observation condition parameter including the measurement data 
and the pupil data; 

second means for calculating an evaluation parameter 

10 indicating quality of visibility by the eye to be examined based 
on the observation condition parameter and/or the obtained 
lower order aberrations and the higher order aberrations; and 
third means for, in accordance with the calculated 
evaluation parameter, (1) estimating one of or two or more of 

15 the visual acuity, the optical characteristic and the 
sensitivity, of the eye to be examined under the observation 
condition corresponding to the environment of a subjective eye 
or ( 2 ) calculating the appropriate correction data suitable for 
the eye to be excunined by changing the lower order aberration. 

20 According to the second solving means of this invention, 

there is provided, an ophthalmic data measurement program for 
causing a computer to execute: 

a first step at which a calculation section obtains 
measurement data indicating a refractive power distribution of 

25 an eye to be examined and pupil data including a value of a pupil 
diameter of the eye to be examined or a pupil diameter image, 
and obtains lower order aberrations and higher order 
aberrations based on an observation condition parameter 
including the measurement data and the pupil data; 

30 a second step at which the calculation section calculates 

5 



an evaluation parameter Indicating quality of visibility by the 
eye to be examined based on the observation condition parameter 
and/ or the obtain lower order aberrations and the higher order 
aberrations ; and 
5 a third step at which in accordance with the calculated 

evaluation pareuneter, the calculation section estimates one of 
or two or more of a visual acuity, an optical characteristic 
and a sensitivity of the eye to be examined under an observation 
condition corresponding to an environment of a subjective eye, 

10 or calculates appropriate correction data suitable for the eye 
to be examined by changing the lower order aberration. 

According to the third solving means of this invention, 
there is provided, an ophthalmic data measurement program for 
causing a computer to execute: 

15 a first step at which a calculation section receives 

measurement data indicating a refractive power distribution of 
an eye to be examined, and obtains lower order aberrations and 
higher order aberrations based on the measurement data; 

a second step at which the calculation section calculates 

20 an evaluation parameter indicating quality of visibility by the 
eye to be examined based on the obtained lower order aberrations 
and the higher order aberrations; and 

a third step at which the calculation section calculates 
appropriate correction data suitable for the eye to be examined 

25 by changing the lower order aberration in accordance with the 
calculated evaluation pareuneter. 

According to the fourth solving meeuis of this invention, 
there is provided, an eye characteristic measuring apparatus 
comprising: 

30 a first light source part to emit a light flux of a first 
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wavelength ; 

a first illtuninating optical system for performing 
illumination to condense the light flux from the first light 
source part on a vicinity of a retina of an eye to be exeunined; 
5 a first light receiving optical system for receiving a 

part of the light flux reflected by and returning from the retina 
of the eye to be examined through a first conversion member to 
convert it into at least substantially 17 beams; 

a first light receiving part for receiving the received 
10 light flux of the first light receiving optical system; and 

a calculation section for receiving pupil data including 
a pupil image of the eye to be examined in a measurement 
environment , calculating a pupil diameter under the measurement 
environment based on the received pupil data, and obtaining an 
15 optical characteristic of the eye to be examined based on the 
calculated pupil diameter and an output signal from the first 
light receiving part. 



Brief Description of the Drawings 
20 Fig, 1 is a structural view of an optical system 100 of 

an eye optical characteristic measuring apparatus . 

Fig. 2 is a structural view of an electric system 200 of 
the eye optical characteristic measuring apparatus. 

Fig. 3 is an explanatory view of a Landolt ' s ring. 
25 Fig. 4 is flowchart of ophthalmic data measurement. 

Fig. 5 is a sub-flowchart concerning calculation of a 
pupil diameter and measurement of eye optical system data. 

Fig. 6 is an explanatory view of pupil diameter 
calculation. 

30 Fig. 7 is a first flowchart of simulation of visual 
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acuity . 

Fig. 8 Is a flowchart of simulation of a retinal Image. 
Fig. 9 Is an explanatory view of template matching. 
Fig. 10 Is a flowchart of Landolt's ring template 
5 matching . 

Pig. 11 Is a second flowchart (1) of simulation of visual 
acuity. 

Fig. 12 is a second flowchart (2) of simulation of visual 
acuity. 

10 Fig. 13 is a third flowchart of simulation of visual 

acuity . 

Fig. 14 is a fourth flowchart of simulation of visual 
acuity. 

Fig. 15 is an explanatory view of contrast sensitivity. 
15 Fig. 16 shows a display excunple of visual acuity 

estimation by template matching. 

Fig. 17 shows a display example concerning comparison 
between pre-compensatlon and post -compensation. 

Fig. 18 is an expleuiatory view of an example of 
20 prescription data for eyeglasses/contacts. 

Fig. 19 is an explanatory view of an example of data for 
refractive surgery. 

Fig. 20 is an explanatory view of an example of 
prescription data for eyeglasses/contacts when an 
25 environmental condition is changed. 

Fig. 21 is an explanatory view of an example of pupil data 
when an environmental condition is changed. 

Fig. 22 is a comparison view of prescription data for 
eyeglasses/contacts with respect to measurement with a constant 
30 pupil diameter. 



Fig, 23 Is a flowchart of correction image simulation. 

Fig. 24 is a flowchart concerning a first example of best 
image condition calculation. 

Fig. 25 is a flowchart concerning a second excunple of best 
5 image condition calculation. 

Fig. 26 is a view showing a display example of best image 
display - Strehl ratio optimization. 

Fig. 27 is a view showing a display example of best image 
display - PTF optimization. 
10 Fig. 28 is a view showing a display example concerning 

comparison between pre -compensation and post -compensation. 

Fig. 29 is an explanatory view of an example of 
prescription data for eyeglasses/contacts. 

Fig. 30 is an explanatory view of an example of data for 
15 refractive surgery. 

Fig. 31 is an explanatory view of an example of 
prescription data for eyeglasses/contacts when an 
environmental condition is changed. 

Fig. 32 is a comparison view of prescription data for 
20 eyeglasses/contacts with respect to measurement with a constant 
pupil diameter. 

Best Mode for Carrying Out the Invention 

25 Fig. 1 is a structural view of an optical system 100 of 

an eye optical characteristic measuring apparatus (an 
ophthalmic data measuring apparatus ) . 

The optical system 100 of the eye optical characteristic 
measuring apparatus is an apparatus for measuring an optical 

30 characteristic of an eye 60 to be measured as an object, and 
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Includes a first illuminating optical system 10, a first light 
receiving optical system 20, a second light receiving optical 
system 30, a common optical system 40, an adjusting optical 
system 50, a second illuminating optical system 70, and a second 
5 light sending optical system 80 • Incidentally, with respect to 
the eye 60 to be measured, a retina 61 and a cornea 62 are shown 
in the drawing. 

The first illiiminating optical system 10 includes, for 
example, a first light source part 11 for emitting a light flux 

10 of a first wavelength, and a condensing lens 12, and is for 
illuminating a minute area on the retina (retina) 61 of the eye 
60 to be measured with the light flux from the first light source 
part 11 so that its illumination condition can be suitably set. 
Incidentally, here, as an exeimple, the first wavelength of the 

15 illuminating light flux emitted from the first light source part 
11 is a wavelength (for example, 780 nm) of an infrared range. 

Besides, it is desirable that the first light source part 
11 has a high spatial coherence and a low temporal coherence. 
Here, the first light source part 11 is, for example, a super 

20 luminescence diode (SLD) , and a point light source having high 
luminescence can be obtained. Incidentally, the first light 
source part 11 is not limited to the SLD, and for example, a 
laser having a high spatial coherence and a high temporal 
coherence can also be used by inserting a rotation diffused 

25 plate or declination prism (D prism) or the like to suitably 
lower the temporal coherence. Further, an LED having a low 
spatial coherence and a low temporal coherence can also be used, 
if light quantity is sufficient, by inserting, for example, a 
pinhole or the like at a position of a light source in an optical 

30 path. 



The first light receiving optical system 20 includes, for 
example, a collimator lens 21, a Hartmann plate 22 as a 
conversion member for converting a part of a light flux (first 
light flux) reflected and returned from the retina 61 of the 
5 eye 60 to be measured into at least 17 beams, and a first light 
receiving part 23 for receiving the plural beeims converted by 
the Hartmann plate 22, and is for guiding the first light flux 
to the first light receiving part 23. Besides, here, a CCD with 
little readout noise is adopted for the first light receiving 

10 part 23, and as the CCD, a suitable type of CCD, for exeimple, 
a general low noise type of CCD, a cooling CCD of 1000 ^ 1000 
elements for measurement, or the like is applicable. 

The second illuminating optical system 70 includes a 
second light source 72 and a Placido's disk 71. Incidentally, 

15 the second light source 72 can be omitted. The Placido's disk 
(PLACIDO'S DISK) 71 is for projecting an index of a pattern 
composed of plural co- axial rings. Incidentally, the index of 
the pattern composed of the plural co- axial rings is an example 
of an index of a specified pattern, and a different suitable 

20 pattern can be used. Then, after an alignment adjustment 
described later is completed, the index of the pattern composed 
of the plural co-axial rings can be projected. 

The second light sending optical system 80 is for mainly 
performing, for example, the alignment adjustment described 

25 later, and measurement and adjustment of a coordinate origin 
and a coordinate eixis, and includes a second light source part 
31, a condensing lens 32, and a beam splitter 33. 

The second light receiving optical system 30 includes a 
condensing lens 34 and a second light receiving part 35. The 

30 second light receiving optical system 30 guides a light flux 
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(second light flux), which is originated from the pattern of 
the Placldo's disk 71 illuminated from the second illuminating 
optical system 70 and is reflected and returned from the 
anterior eye part or the cornea 62 of the eye 60 to be measured, 
5 to the second light receiving part 35 • Besides , it can also guide 
a light flux, which is emitted from the second light source part 
31 and is reflected and returned from the cornea 62 of the eye 
60 to be measured, to the second light receiving part 35. 
Incidentally, as the second wavelength of the light flux emitted 

10 from the second light source part 31, for exeimple, a wavelength 
different from the first wavelength (here, 780 nm) and longer 
(for example, 940 nm) than that can be selected. 

The common optical system 40 is disposed on an optical 
axis of the light flux emitted from the first illuminating 

15 optical system 10, can be included in the first and the second 
illximinating optical systems 10 and 70, the first and the second 
light receiving optical systems 20 and 30, the second light 
sending optical system 80 and the like in common, and includes, 
for example, an afocal lens 42, beam splitters 43 and 45, and 

20 a condensing lens 44. The beam splitter 43 is formed of such 
a mirror (for example, a dichroic mirror) that the wavelength 
of the second light source part 31 is sent (reflected) to the 
eye 60 to be measured, and the second light flux reflected and 
returned from the retina 61 of the eye 60 to be measured is 

25 reflected, and on the other hand, the wavelength of the first 
light source part 11 is transmitted. The beam splitter 45 is 
formed of such a mirror (for example, a polarization beam 
splitter) that the light flux of the first light source part 
11 is sent (reflected) to the eye 60 to be measured, and the 

30 first light flux reflected and returned from the retina 61 of 



the eye 60 to be measured Is transmitted • By the becun splitters 
43 and 45, the first and the second light fluxes do not mutually 
enter the other optical systems to generate noise. 

The adjusting optical system 50 Is for mainly performing, 
5 for example, a working distance adjustment described later. 
Includes a third light source part 51, a fourth light source 
part 55, condensing lenses 52 and 53, and a third light receiving 
part 54, and Is for mainly performing the working distance 
adjustment . 

10 A third Illuminating optical system 90 Includes an 

optical path for projection of an Index for causing, for excimple, 
fixation of the eye 60 to be measured or fogging, and Includes 
a fifth light source part (for example, a lamp) 91, a fixation 
target 92 and a relay lens 93. The fixation target 92 can be 

15 Irradiated to the retina 61 by the light flux from the fifth 
light source part 91, and the eye 60 to be measured Is made to 
observe Its Image. The fixation target 92 and the retina 61 are 
put In a conjugated relation by the third Illuminating optical 
system 90. Besides, the fifth light source part 91 Is also a 

20 light source (anterior ocular segment Illuminating part) to 
Illuminate the anterior ocular segment of the eye 60 to be 
measured under different luminance conditions. The light 
amount of the fifth light source part 91 Is adjusted, so that 
the Illumination state of the eye 60 to be measured Is changed 

25 and the size of the pupil can be changed. Incidentally, as the 
anterior ocular segment Illuminating part. In addition to the 
fifth light source part 91, an appropriate light source such 
as the second light source 72 may be used. 

Next, the alignment adjustment will be described. The 

30 alignment adjustment Is mainly carried out by the second light 
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receiving optical system 30 and the second light sending optical 
system 80. 

First, the light flux from the second light source part 
31 illuminates the eye 60 to be measured as the object with the 
5 substantially parallel light flux through the condensing lens 
32, the beam splitters 33 and 43, and the afocal lens 42. The 
reflected light flux reflected by the cornea 62 of the eye 60 
to be measured is emitted as a divergent light flux such as is 
emitted from a point at the half of the radius of curvature of 
10 the cornea 62. The divergence light flux is received as a spot 
image by the second light receiving part 35 through the afocal 
lens 42, the beam splitters 43 and 33, and the condensing lens 
34. 

Here , in the case where the spot image on the second light 

15 receiving part 35 is outside the optical axis, the main body 
of the eye optical characteristic measuring apparatus is moved 
and adjusted vertically and horizontally, and the spot image 
is made to coincide with the optical axis. As stated above, when 
the spot image coincides with the optical axis, the alignment 

20 adjustment is completed. Incidentally, with respect to the 
alignment adjustment, the cornea 62 of the eye 60 to be measured 
is Illuminated by the third light source part 51, and an image 
of the eye 60 to be measured obtained by this illumination is 
formed on the second light receiving part 35, and accordingly, 

25 this image may be used to make the pupil center coincide with 
the optical axis. 

Next, the working distance adjustment will be described. 
The working distance adjustment is mainly carried out by the 
adjusting optical system 50. 

30 First, the working distance adjustment is carried out by, 
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for example. Irradiating the eye 60 to be measured with a 
parallel light flux emitted from the fourth light source part 
55 and close to the optical axis, and by receiving the light 
reflected from the eye 60 to be measured through the condensing 
5 lenses 52 and 53 by the third light receiving part 54. Besides, 
in the case where the eye 60 to be measured is in a suitable 
working distance, a spot image from the fourth light source part 
55 is formed on the optical eixis of the third light receiving 
part 54. On the other hand, in the case where the eye 60 to be 

10 measured goes out of the suitable working distance, the spot 
image from the fourth light source part 55 is formed above or 
below the optical axis of the third light receiving part 54* 
Incidentally, since the third light receiving part 54 has only 
to be capable of detecting a change of a light flux position 

15 on the plane containing the fourth light source part 55, the 
optical axis and the third light receiving part 54, for example, 
a one -dimensional CCD arranged on this plane, a position sensing 
device (PSD) or the like is applicable. 

Next, a positional relation between the first 

20 illuminating optical system 10 and the first light receiving 
optical system 20 will be described. 

The beam splitter 45 is inserted in the first light 
receiving optical system 20, and by this beam splitter 45, the 
light from the first illuminating optical system 10 is sent to 

25 the eye 60 to be measured, and the reflected light from the eye 
60 to be measured is transmitted. The first light receiving part 
23 included in the first light receiving optical system 20 
receives the light transmitted through the Hartmann plate 22 
as the conversion member and generates a received light signal. 

30 Besides, the first light source part 11 and the retina 
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61 of the eye 60 to be measured form a conjugated relation. The 
retina 61 of the eye 60 to be measured and the first light 
receiving part 23 are conjugate. Besides, the Hartmann plate 
22 and the pupil of the eye 60 to be measured form a conjugated 
5 relation. Further, the first light receiving optical system 20 
forms a substantially conjugated relation with respect to the 
cornea 62 as the anterior eye part of the eye 60 to be measured, 
the pupil, and the Hartmann plate 22. That is, the front focal 
point of the afocal lens 42 is substantially coincident with 

10 the cornea 62 as the anterior eye part of the eye 60 to be measured 
and the pupil. 

Besides, the first illuminating optical system 10 and the 
first light receiving optical system 20 are moved together so 
that a signal peak according to the reflected light at the light 

15 receiving part 23 becomes maximum on the condition that the 
light flux from the first light source part 11 is reflected at 
a point on which it is condensed. Specifically, the first 
illuminating optical system 10 and the first light receiving 
optical system 20 are moved in a direction in which the signal 

20 peak at the first light receiving part 23 becomes large, and 
are stopped at a position where the signal peak becomes maximum. 
By this, the light flux from the first light source part 11 is 
condensed on the eye 60 to be measured. 

Besides, the lens 12 converts a diffused light of the 

25 light source 11 into a parallel light. A diaphragm 14 is 
positioned at an optically conjugated position with respect to 
the pupil of the eye or the Hartmann plate 22. The diaphragm 
14 has a diameter smaller than an effective range of the Hartmann 
plate 22, and the so-called single path aberration measurement 

30 (method in which aberrations of an eye have an influence on only 



the light receiving side) is established. In order to satisfy 
the above, the lens 13 is disposed such that the retina 
conjugated point of the real light beam coincides with the front 
focal position, and further, in order to satisfy the conjugated 
5 relation between the lens and the pupil of the eye, it is disposed 
such that the rear focal position coincides with the diaphragm 
14. 

Besides, after a light beam 15 comes to have a light path 
common to a light beam 24 by the beam splitter 45, it travels 

10 in the same way as the light beam 24 paraxially. However, in 
the single path measurement, the diameters of the light beams 
are different from each other, and the beeun dicuneter of the light 
beam 15 is set to be rather small as compeured with the light 
beam 24. Specifically, the beam diameter of the light beam 15 

15 is, for example, about 1 mm at the pupil position of the eye, 
and the beemi diameter of the light beam 24 can be about 7 mm 
(incidentally, in the drawing, the light beam 15 from the beam 
splitter 45 to the retina 61 is omitted) • 

Next, the Hartmann plate 22 as the conversion member will 

20 be described. 

The Hartmann plate 22 included in the first light 
receiving optical system 20 is a wavefront conversion member 
for converting a reflected light flux into plural beams. Here, 
plural micro-Fresnel lenses disposed on a plane orthogonal to 

25 the optical axis apply in the Hartmann plate 22. Besides, in 
general, with respect to the measurement object part (the eye 
60 to be measured) , in order to measure a sphere of the eye 60 
to be measured, third- order astigmatism aberrations , and other 
higher order aberrations, it is necessary to perform the 

30 measurement with at least 17 beams through the eye 60 to be 



measured . 

Besides, the micro-Fresnel lens Is an optical element, 
and Includes, for example, a ring with a height pitch for each 
wavelength, and a blade optimized for emission parallel to a 
5 condensing point. The micro-Fresnel lens here is subjected to, 
for example, 8 -level optical path length variation employing 
a semiconductor fine working technique, and achieves a high 
condensing efficiency (for example, 98 %) . 

Besides, the reflected light from the retina 61 of the 

10 eye 60 to be measured passes through the afocal lens 42 and the 
collimate lens 21 and is condensed on the first light receiving 
part 23 through the Hartmann plate 22 . Accordingly, the Hartmann 
plate 22 includes a wavef ront conversion member for converting 
the reflected light flux into at least 17 beams. 

15 Fig. 2 is a structural view of an electrical system 200 

of the eye optical characteristic measuring apparatus . The 
electrical system 200 of the eye optical characteristic 
measuring apparatus includes, for example, a calculation 
section 210, a control section 220, a display section 230, a 

20 memory 240, an input section 270, a first driving section 250, 
a second driving section 260, and a third driving section 280. 
Besides, the calculation section 210 can include a pupil data 
formation section 215, an image data formation section 211, a 
judgment section 212, and a correction element specifying 

25 section 123. 

The calculation section 210 receives a received light 
signal (4) obtained from the first light receiving part 23, a 
received light signal (7) obtained from the second light 
receiving part 35, and a received light signal (10) obtained 

30 from the third light receiving part 54, and performs an 



arithmetical operation on the origin of coordinates, coordinate 
eixis, movement of coordinates, rotation, pupil diameter , ocular 
aberrations, corneal wavefront aberrations, Zernike 
coefficients, aberration coefficients, visual acuity 
5 simulation, Strehl ratio (Strehl ratio), phase shift (PTF, 
phase shift) , white light MTF, Landolt's ring pattern, contrast 
sensitivity and the like. The processing of ideal observed 
analysis as proposed by Wilson Geiseler may be performed 
(Geisler, W.S. 1989 Psychological Review 96, pp. 267-324). 

10 Besides, signals corresponding to such calculation results are 
output ted to the control section 220 for performing the whole 
control of an electric driving system, the display section 230 
and the memory 240, respectively. Incidentally, the details of 
the calculation section 210 will be described later. 

15 The pupil data formation section 215 forms pupil data from 

an anterior ocular segment image. For example, the pupil data 
formation section 215 receives the anterior ocular segment 
image from the second light receiving part 35, calculates points 
on an edge of a pupil, a focal point, a major axis and a minor 

20 axis when the pupil is elliptic, and obtains a pupil diameter. 
When the pupil area shape is not circular but elliptic or is 
still another shape, this is specified to obtain measured values 
used for analysis. 

In correction data calculation in after-mentioned 

25 template matching or visual acuity simulation, based on the 
measurement data indicating at least wavefront aberrations of 
the eye to be examined and in view of a correction element for 
refraction correction, the image data formation section 211 
performs the simulation of visibility of the index for eye 

30 examination, and forms index retinal image data. The wavefront 
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aberrations of the eye to be examined include higher order 
aberrations. That is, parameters of distributions concerning 
all refractions are included. Based on the index retinal image 
data formed by the image data formation section 211, the 
5 Judgment section 212 judges whether or not the index for eye 
examination is seen. 

Besides, the correction element specifying section 213 
specifies the correction element to be given to the image data 
formation section 211. Further, based on the index retinal 

10 image data which is corrected with the correction element 
specified by the correction element specifying section 213 and 
is formed by the image data formation section 211, the judgment 
section 212 judges whether or not the appropriate correction 
element is specified. Besides, the correction element 

15 specifying section 213 specifies a correction element based on 
the result of the judgment section 212, and repeatedly changes 
the correction element until the judgment section 212 judges 
that it is the appropriate correction element. The correction 
element is one of or a combination of two or more of a spherical 

20 power, an astigmatic power, and an astigmatic axis angle. 

The control part 220 controls lighting and extinction of 
the first light source part 11 on the basis of the control signal 
from the arithmetic part 210, or controls the first driving part 
250 and the second driving part 260. For example, on the basis 

25 of the signals corresponding to the operation results in the 
arithmetic part 210, the control part outputs a signal (1) to 
the first light source part 11, outputs a signal (5) to the 
Placido's disk 71, outputs a signal (6) to the second light 
source part 31, outputs a signal (8) to the third light source 

30 part 51, outputs a signal (9) to the fourth light source part 
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55, outputs a signal (11) to the fifth light source part 91, 
and outputs signals to the first driving part 250 and the second 
driving part 260. 

The first driving part 250 is for moving the whole first 
5 illuminating optical system 10 in the optical axis direction 
on the basis of, for excunple, the received light signal (4) 
inputted to the arithmetic part 210 from the first light 
receiving part 23, and outputs a signal (2) to a not-shown 
suitable lens movement means and drives the lens movement means . 

10 By this, the first driving part 250 can perform the movement 
and adjustment of the first illuminating optical system 10, 
The second driving part 260 is for moving the whole first 
light receiving optical system 20 in the optical axis direction 
on the basis of, for excunple, the received light signal (4) 

15 inputted to the arithmetic part 210 from the first light 
receiving part 23, and outputs a signal (3) to a not -shown 
suitable lens movement means , and drives the lens movement means . 
By this, the second driving part 260 can perform the movement 
and adjustment of the first light receiving optical system 20. 

20 The third driving section 280 is for moving a fixation 

index 92 of the third illuminating optical system 90, and 
outputs a signal (12) to a not-shown appropriate movement means 
and drives the movement means. By this, the third driving 
section 280 can perform the movement euid adjustment of the 

25 fixation index 92 of the third illuminating optical system 90. 

2. Zernike analysis 

Next, a Zernike analysis will be described. A generally 
30 known method of calculating Zernike coefficients Ci^^'^ from 
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Zernike polynomials will be described. The Zernike coefficients 
Ci^^'^ are important parameters for grasping the optical 
characteristic of the subject eye 60 on the basis of inclination 
angles of the light f Ixixes obtained by the first light receiving 
part 23 through the Hartmann plate 22. 

Wavefront aberrations W(X, Y) of the subject eye 60 are 
expressed using the Zernike coefficients Ci^^"^ and the Zernike 
polynomials Z±^^'^ by the following expression. 

1=0 y=o 

(1) 

Where, (X, Y) denotes vertical and horizontal coordinates 
of the Hartmann plate 22. 

Besides, with respect to the wavefront aberrations W(X, 
Y) , when the horizontal and vertical coordinates of the first 
light receiving part 23 are denoted by (x, y) , a distance between 
the Hartmann plate 22 and the first light receiving part 23 is 
denoted by f , and a movement distance of a point image received 
by the first light receiving part 23 is denoted by (Ax, Ay), 
the following expression is established. 

dW(X,Y) ^Ax 
dX " / 

(2) 

^W(XJ) ^^y 

dY ^ f 

(3) 

Where, the Zernike polynomials Zi^^"^ are expressed by the 
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5 



10 



following nvunerlcal expressions ( 4 ) and ( 5 ) . (More specifically 
expressions, for exeunple, see JP-A-2002-209854 . ) 



[cos J 



m>0 sin 

m^O cos 



(4) 



in-m)/2 r»\ i 

Kir)- Z (-0'— r ^^ji^ r'-' 



(5) 

Incidentally, with respect to the Zernlke coefficients 
Ci^^'^, specific values can be obtained by minimizing the squared 
error expressed by the following numerical expression. 



data number 



\ dX f\ \ dY f\ 



(6) 

Where, W(X, Y) : wavefront aberrations, (X, Y) : Hartmann 
plate coordinates, (Ax, Ay): a movement distance of a point 
Image received by the first light receiving part 23, f: a 
15 distance between the Hartmann plate 22 and the first light 
receiving part 23. 

The arithmetic part 210 calculates the Zernlke 
coefficients Ci^^"^, and uses this to obtain eye optical 
characteristics such as spherical aberrations, coma 
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aberrations, and astigmatism aberrations. 



(Normalization of Pupil Diameter) 

A Zernike polynomial always indicates a shape in a circle 
5 with a radius of 1, and when Zernike analysis is performed at 
a certain pupil diameter (pupil diameter) , normalization is 
performed with the pupil radius. For example, when the center 
coordinate of the pupil with the pupil radius rp is made (0, 
0) , a point P(X, Y) in the pupil is made P(X/rp, Y/rp) when the 

10 Zernike analysis is performed. When the barycentric point of 
a spot of a Hartmann image is P, a reference lattice point Pref 
(Xref# Yref) Corresponding to this point is made Pre£(Xref /^p, 
Yref/i^p) # and a movement distance of a point image is obtained, 
and the Zernike coefficients are calculated. The actual 

15 wavefront (wavefront in which the coordinates are not 
normalized) W(X, Y) is expressed by a following expression. 

1=0 7=0 /=o y=o 

(7) 

Where, (X, Y) : coordinates not normalized, (Xg, ys) : coordinates 
20 normalized. 



3. Landolt*s ring 



Fig. 3 is an explanatory view of a Landolt*s ring. 
25 Hereinafter, preparation of data of a luminous 

distribution function Land(x, y) of the Landolt*s ring will be 
described. Fig. 3 shows Landolt's ring of high contrast in upper 
stand, and Landolt's ring of low contrast in lower stand. 
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The Landolt ' s ring is expressed by the reciprocal of a 
recognizable minimvun visual angle, and the ability to be capable 
of recognizing a visual angle of one minute is called visual 
acuity of 20/20 . For example, if the recognizable minimum visual 
5 angle is 2 minutes, the visual acuity is defined as 20/40, and 
if 10 minutes , the visual acuity is defined as 20/200 . In general, 
the Landolt ' s ring uses , as an index, a ring in which a gap being 
1/5 of the size of the outside ring is provided as shown in the 
drawing . 

10 When the visual acuity is V, the size d of the Landolt 's 

ring projected on the retina is calculated by 



t/=5x2/?tan 



(8) 

(R: a distance between a pupil and an image point (retina)) 
15 On the basis of this expression and the definition of the 

Landolt 's ring, a black portion of the Landolt 's ring is made 
0 (or 1), a white portion thereof is made 1 (or 0), and the 
luminous distribution function Land(x, y) of the Landolt 's ring 
is prepared. The data of the prepared luminous distribution 
20 function Land(x, y) is stored in the memory 240, is read out 
by the arithmetic part 210, and is set correspondingly to 
predetermined visual acuity. 

A high contrast original image is such that for example, 
the contrast of the black portion and the white portion of a 
25 Landolt 's ring is 100% (for example, the white is 0 and the black 
is 1), or the figure portion of a Landolt 's ring is black (10 
cd/m^ or less) and the background is white (100 cd/m^) and the 
actual contrast is 90% or more. Here, the contrast is such that 
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Mlchelson contrast (I white - I black) /(I white + I black) is 
expressed by % . On the other hand, as a low contrast original 
image, one in which the contrast of the black portion and the 
white portion of a Landolt ' s ring is 10% (for example, the white 
5 is 0 and the black is 0.1) can be used. These contrasts have 
an accuracy of approximately ±1%. Incidentally, in addition 
to this example, an appropriate contrast original image may be 
used. As a luminance distribution function Land (x, y) stored 
in the memory 240, a high contrast one and a low contrast one 
10 are respectively formed and are stored. 

4. Ophthalmic data measuring method 

Fig. 4 is a flowchart of ophthalmic data measurement. 
First, the eye optical characteristic measuring 

15 apparatus makes alignment of X, Y and Z axes of the pupil position 
of the eye 60 to be measured (SlOl). Next, the measuring 
apparatus moves the origin of a movable section (S103). For 
example, the Hartmann plate 22, the Placido*s disk 71 or the 
like is set to zero diopter. The calculation section 210 

20 measures the data of the eye optical system, such as the pupil 
diameter, ocular aberrations and Zernike coefficients, on the 
basis of the measured received light signals (4), (7) and/or 
(10) (S105). The calculation section 210 performs a visual 
acuity simulation or a correction image simulation (S107). 

25 In the visual acuity simulation, for example, at step S107, 

the calculation section 210 uses, as an evaluation parameter 
indicating the quality of visibility by the eye 60 to be examined, 
a comparison result between a simulation result of visibility 
of the index for eye examination and a predetermined template 

30 and/or an MTF (Modulation Transfer Function) indicating the 
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transfer characteristic of the eye to be examined, and estimates 
the visual acuity of the eye to be examined or the sensitivity 
in accordance with the evaluation parameter. Incidentally, as 
the visual acuity^ the index for eye examination is suitably 
5 set, so that the high contrast visual acuity and low contrast 
visual acuity can be estimated. Besides, the calculation 
section 210 estimates optical characteristics such as the MTF 
of the eye to be examined and the point spread function (PSF) . 
Besides, in the correction image simulation, for example, 

10 the calculation section 210 obtains appropriate correction data 
while using one of or two or more of the Strehl ratio, the PTF, 
and the MTF (Modulation Transfer Function) as the evaluation 
parameter indicating the quality of visibility by the eye 60 
to be measured. Besides , the calculation section 210 may obtain 

15 appropriate correction data by, for exeimple, performing the 
simulation of visibility of the index for eye examination and 
using the comparison result to a predetermined template as an 
evaluation pareimeter. 

Incidentally, the details of step S105 and S107 will be 

20 described later. The calculation section 210 outputs data to 
the display section 230 and the memory 240 ( S109 ) . Incidentally, 
in the case where data output has already been made in the former 
processing, the processing of step S109 may be omitted. 

Fig. 5 is a sub-flowchart concerning calculation of the 

25 pupil diameter at step S105 and measurement of data of the eye 
optical system. Besides, Fig. 6 is an explanatory view of pupil 
diameter calculation. 

First, the calculation section 210 acquires a Hartmann 
image and an anterior ocular segment image from the first light 

30 receiving part 20 and the second light receiving part 35 (S601) . 
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The calculation section 210 causes the fifth light source part 
91 to Illuminate the eye 60 to be measured In an Illumination 
state of a desired environmental condition (observation 
condition) , and acquires the Hartmann Image and the anterior 
5 ocular segment Image from the first light receiving part 20 and 
the second light receiving part 35. For example, the 
calculation section 210 causes the display section 230 to 
display Instructions to select an environmental condition under 
which the visual acuity or sensitivity Is estimated, and the 

10 selected environmental condition may be Inputted from the Input 
section 270. The environmental condition Includes, for 
example, "seeing in the daytime", "seeing in the twilight" 
"seeing in a room (under a fluorescent lamp)", "seeing in the 
nighttime", "normal visual acuity measurement" and the like. 

15 Next, the calculation section 210 refers to, for example, a 
table which is previously stored in the memory 240 and in which 
the environmental conditions and the illumination states 
correspond to each other, and acquires the illumination state 
corresponding to the inputted environmental condition. The 

20 illumination states under the respective environmental 
conditions can be made such that for example, the case of "normal 
visual acuity measurement" is 50 [Ix] , "seeing in the daytime" 
is 100000 [Ix] , and "in a room (fluorescent lamp) " is 2000 [Ix] . 
Incidentally, with respect to these values, em appropriate 

25 value corresponding to the environmental condition can be used. 
As the environment, it is desirable to use a fixation target 
larger than a normal one. Here, although the eye 60 to be 
examined is illuminated in the illumination state of the desired 
environmental condition by the fifth light source part 91, a 

30 structure may be made such that the illumination state is formed 
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by using the surrounding illumination of the eye to be examined 
or the background illumination. 

The calculation section 210 outputs a signal (11) 
corresponding to the acquired illumination state to the fifth 
5 light source part 91 through the control section 220, and causes 
the eye 60 to be measured to be illuminated. Besides, the 
calculation section 210 sequentially changes the illumination 
state from a dark one to a bright one, and can acquire the 
Hartmann images and anterior ocular segment images in the plural 

10 illumination states. 

Incidentally, the calculation section 210 may omit step 
S601, and reads Hartmann image data previously measured and 
stored in the memory 240, and pupil data including one of the 
anterior ocular segment image, the pupil shape such as points 

15 on the pupil edge, and pupil diameter. Besides, for example, 
the calculation section 210 may acquire the anterior ocular 
segment image by reading, as the pupil data in an electric carte, 
photographic data photographed in the past and stored in the 
memory 240 from the memory. 

20 Next, based on the acquired anterior ocular segment image, 

the calculation section 210 detects, for example, 36 (n - 36) 
points Pi (i = 1 to n) on the edge of the pupil (S603). The 
calculation section 210 detects the change (light and shade on 
the image) of the acquired light amount of the anterior ocular 

25 segment image by a method of image processing, and can obtain 
points on the edge of the pupil. In Fig. 6, the detection points 
Pi are points indicated by marks of "+". 

Next, the calculation section 210 performs elliptic 
fitting which is fittest to the detected points on the edge of 

30 the pupil (S605). First, the calculation section 210 obtains 



the focal points (points Fl and F2 in Fig. 6) of the ellipse • 
For example, the calculation section 210 reads the coordinates 
of two points previously set as the initial values of the focal 
points from the memory 240. Next, the calculation section 210 
obtains distances from the detection point Pi to the two read 
points , and the sum of the distances is made Li. The calculation 
section 210 obtains the sum Li of the distances concerning all 
the detection points Pi, and obtains a mean value A of Li. 
Further, the calculation section 210 uses a method of the least 
square approximation or the like to calculate two points where 
a square error Se of the sum Li of the distance and the mean 
value A expressed by a following expression becomes minimum, 
and consequently, the focal points of the ellipse can be 
obtained . 



Where, Li: the sum of distances from the point Pi on the edge 
to the two points Fl and F2, A: the mean value of Li at the 
respective points on the edge, n: the number of detected points 
on the edge. Incidentally, the focal points of the ellipse may 
be obtained by an appropriate method other than this . 

Next, the calculation section 210 obtains the sum L of 
distances from one point on the ellipse to the focal points. 
Incidentally, the calculation section 210 may make the 
foregoing mean value A the sum L of the distances from one point 
on the ellipse to the focal points. Next, the calculation 
section 210 calculates the pupil diameter from the length (major 
axis) of the long axis of the ellipse and the length (minor axis) 



n 




(9) 



/=1 
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of the short clxIs (S607). The length 2a of the long axis and 
the length 2b of the short axis can be expressed by following 
expressions • 

2a = L 
26 = 2. 



2) 



-c 



= 2. 



4 4 



(10) 



Where, L: the sum of distances from a point on the edge to the 
focal points , (xl , yl ) , (x2 , y2 ) : the focal points of the ellipse . 
When it is assvmied that the pupil diameter dp is, for exeunple, 
the mean value of the length 2a of the long axis and the length 
2b of the short axis, it is expressed by a following expression. 



dp=a + b 



(11) 



Incidentally, instead of making the mean value the pupil 
diameter, an appropriate value based on the length 2a of the 
long axis and the length 2b of the short axis , such as the length 
of the short axis , the length of the long axis , or an intermediate 
15 value of the lengths of the short axis and the long axis, may 
be used. 

The calculation section 210 obtains the pupil center 
position based on, for example, the focal points of the ellipse 
and/or the lengths of the long axis and the short axis, and 
20 further obtains or specifies the limbus center, and may 
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calculate the shift amount of the pupil center position such 
as the shift amount from the limbus center. Besides, the 
calculation section 210 makes the calculated shift amount 
correspond to the pupil diameter and stores it into the memory 
5 240. 

Incidentally, the calculation section 210 may adjust the 
brightness of the fifth light source part 91 to produce the 
illumination state in which the pupil diameter in the 
environment desired by the subjective eye (for example, in an 

10 office, in a classroom, driving in the night, etc. ) is obtained, 
in addition to the illumination state in which the pupil 
diameter in the daytime is obtained. Besides, pupil diameters 
in the above environments are previously measured and may be 
used to perform the analysis . By this , the optimum prescription 

15 value in the environment desired by the subjective eye can be 
analyzed. Incidentally, the calculation section 210 may read 
the measured data and the pupil diameter previously stored in 
the memory 240 instead of the processing of steps S601 to S607. 
The calculation section 210 calculates eye optical system 

20 data based on the pupil dicuneter and the Hartmann image (S609) . 
First, the calculation section 210 detects barycenter points 
of the respective spots from the Hartmann image acquired at step 
S601. Next, the calculation section 210 normalizes the 
barycenter point coordinates detected when the pupil center is 

25 made the origin by the pupil radius rp. Here, the pupil radius 
rp = pupil diameter dp/2. That is, the calculation section 210 
sets the barycenter point Ps(X, Y) in the range of the pupil 
diameter to Ps(X/rp, Y/rp) , and when the barycenter point of the 
spot of the Hartmann image is Ps, the reference lattice point 

30 Pref (Xref . Yref ) Corresponding to this point is made Pref (Xref /rp. 
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Yref/Tp) . The actual wavefront (wavefront in which the 
coordinates are not normalized) W(X, Y) is expressed by a 
following expression. 

/=0 J=0 '=0 7=0 

(12) 

Where, (X, Y) : coordinates not normalized, (Xs, Ys) : coordinates 
normalized. 

The calculation section 210 uses the normalized 
coordinates, and calculates the eye optical system data such 
as the Zernike coefficients and ocular aberrations. Besides, 
the calculation section 210 stores data into the memory 240 at 
an appropriate timing. 

5. Visual acuity simulation 
15 5-1. First flowchart (spherical power correction) of visual 
acuity simulation 

Fig. 7 shows a first flowchart of visual acuity simulation. 
Fig. 7 shows the flowchart in which retinal image simulation 
is performed, a correction spherical power is obtained so that 
20 a Landolt's ring can be detected, and the visual acuity at the 
time of the correction is estimated. Incidentally, in the 
following respective flowcharts, at steps denoted by the same 
reference characters, the same processing is performed. 

First, the calculation section 210 calculates a tentative 
25 spherical power Sr (S1401). As the tentative spherical power 
Sr, for example, a refractive value or a value calculated from 
wavefront aberrations may be used, or a value previously stored 
in the memory 240 or a value inputted from the input section 
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270 may be used. 

Next, the calculation section 210 specifies a spherical 
power Ss for simulation (S1451). Normally, Ss is specified to 
be a weak correction with respect to Sr (for example, Ss = Sr 
5 + 5D) • The calculation section 210 specifies the Landolt ' s ring 
of a previously determined visual acuity Vs (for example, Vs 
=0.1) (S1453). At this time, first, the calculation section 
210 specifies which of high contrast visual acuity and low 
contrast visual acuity is estimated. For example, the 

10 calculation section 210 may specify the high contrast or low 
contrast in accordance with the input from the input section 
270 or the specification previously stored in the memory 240. 
In accordance with the specification, the calculation section 
210 specifies the Landolt 's ring of the high contrast or low 

15 contrast corresponding to the previously determined visual 
acuity Vs. 

The image data formation section 211 of the calculation 
section 210 performs the simulation of a Landolt 's ring retinal 
image to obtain index image data (S1405) . Here, the image data 

20 formation section 211 performs it with respect to the Landolt *s 
ring in a previously determined direction (for example, a gap 
of the ring is provided in an upper, lower, right or left 
direction). That is, in accordance with the wavefront 
aberrations measured at step S105, the image data formation 

25 section 211 obtains index image data indicating the visibility 
of the Landolt 's ring by simulation. The specific processing 
of this simulation will be described later. 

Next, the judgment section 212 of the calculation section 
210 performs Landolt 's ring template matching (S1407). The 

30 Judgment section 212 performs the template matching between the 
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Index image data obtained by the simulation and the Landolt ' s 
ring in a certain direction, and stores the direction at that 
time and the score n indicating a coincidence degree into the 
memory 240. The specific processing thereof will be described 
5 later . 

The judgment section 212 judges whether the template 
matching is performed in all directions (S1409) • Here, in the 
case of No, advance is made to step S1407, and the processing 
is repeated until the template matching is performed in all the 

10 directions • On the other hand, in the case of Yes at step S1409 , 
the judgment section 212 judges whether the highest score nh 
of the score n matches the direction of the Landolt 's ring of 
the index image data simulated at step S1405 (S1411). Here, 
in the case of Yes, the judgment section 212 judges whether the 

15 score nh is higher than a previously determined threshold value 
in the memory 240 or the like (S1413). Incidentally, as the 
threshold value (threshold value by which a judgment is made 
as to whether the Landolt 's ring can be discriminated), for 
example, a value obtained by comparison with subjective values 

20 of many normal eyes in the past can be used. 

In the case of No at step S1411 or S1413, the judgment 
section 212 judges whether Ss exceeds a previously determined 
allowable value (for example, Sr-5D) (S1415). Here, in the case 
of No, the correction element specifying section 213 sets the 

25 correction element of Ss to be slightly strong (for example, 
SS-0.25D) (S1417), and the image data formation section 211 
performs the simulation of the Landolt 's ring retinal image 
based on this correction element. The calculation section 210 
performs the processing subsequent to step S1407 concerning the 

30 index image data obtained by this simulation. On the other hand. 



In the case of Yes at step S1415, the judgment section 212 judges 
that the Landolt's ring can not be detected (S1419) , and stores 
the direction at this time and that the detection could not be 
made in this direction into the memory 240. 
5 After step S1419 or in the case of Yes at step S1413, the 

judgment section 212 judges whether simulation has been 
performed in all directions of the Landolt ' s ring ( S1421 ) • Here , 
in the case of No, return is made to step S1405, and the 
calculation section 210 repeats the foregoing processing in all 

10 directions . On the other hand, in the case of Yes at step S1421 , 
the judgment section 212 judges whether detection could be made 
in number equal to or more than half of specified number of 
directions ( S1455 ) • 

In the case of Yes at step S1455, the correction element 

15 specifying section 213 specifies S = Ss and V = Vs , and specifies 
the Landolt 's ring of visual acuity Vs = Vs + 0.1 (S1457). At 
this time, as the specified Landolt 's ring, in accordance with 
the specification at step S1453, the Landolt 's ring of the high 
contrast or low contrast is specified. Thereafter, advance is 

20 made to step S1405, the image data formation section 211 
performs the simulation of the retinal image based on the 
specified correction element and Landolt ' s ring to obtain index 
image data, and performs the processing subsequent to step S1407. 
On the other hand, in the case of No at step S1455, the 

25 calculation section 210 performs data output (S1423) . That is, 
the calculation section 210 displays the visual acuity V at this 
time , the spherical power S = Ss , the direction of the Landolt ' s 
ring which could be detected, simulation result and the like 
on the display section 230 and stores them into the memory 240. 

30 Incidentally, the calculation section 210 may use decimal 



visual acuity as the visual acuity, or may use logMAR (log 
Minimum Angle Resolution) visual acuity. The logMAR visual 
acuity is the visual acuity expressed by the logarithm of 
minimum vision. Incidentally, the data to be displayed and 
5 stored are not limited to the foregoing, and appropriate data 
can be displayed and stored. Besides, it may be suitably 
selected among the foregoing data. For exeunple, data except 
the visual acuity V may be displayed. 



10 at the above step S1405. First, the calculation section 210 
calculates a pupil function f(x, y) based on the wavefront 
aberration W(X, Y) obtained at step S105 of Fig. 4 and the 
specified correction element by a following expression (S204) • 



calculates the luminance distribution function Land(x, y) of 
the Landolt's ring (or an arbitrary image) (S205). The 
calculation section 210 performs the two -dimensional Fourier 
transform of Land(x, y) to obtain a spatial frequency 

20 distribution FR(u, v) (S207). The calculation section 210 
calculates a spatial frequency distribution OTF of the eye based 
on the pupil function, and obtains a frequency distribution OR(u, 
v) after passing through the eye optical system by multiplying 
the spatial frequency distribution FR(u, v) of the Landolt's 

25 ring (or arbitrary image) by the spatial frequency distribution 
OTF(u, v) of the eye (S209) as follows. 

FR(u, v) X OTF(u, v) ^ OR(u, v) 
Incidentally, the specific calculation method of the OTF will 
be described later. 



Fig. 8 is a flowchart of the simulation of a retinal image 




(13) 



15 



The calculation section 210 refers to the memory 240 and 
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Next, the calculation section 210 performs the 
two-dimensional inverse Fourier transform of OR(u, v) to obtain 
the luminance distribution image Landlmage (X, Y) of the 
Landolt's ring (or arbitrary image) (S211). 
5 Fig. 9 is an explanatory view of template matching of step 

S1407. As shown in the drawing, a template image (lower 
drawing) is specified correspondingly to the Landolt's ring 
original image (upper drawing) , and the template image as stated 
above is stored in the memory 240 correspondingly to an 

10 identifier indicating the size of the Landolt ' s ring. Although 
the template image is such that in this example, b is made b 
- 1.5a, the number of pixels of the Landolt 's ring is Nl, the 
pixel number is 1 , the number of pixels of a blurred point image 
part around the Landolt 's ring is N2, and the pixel value is 

15 -N1/N2, it is not limited to this but can be suitably specified. 
Besides, although the Landolt "s ring original image shown on 
the upper part of Fig • 9 shows the Landolt ' s ring original image 
of high contrast, also in the case where the Landolt 's ring 
original image of low contrast is used, a similar template can 

20 be used. 

Fig. 10 shows a flowchart of the Landolt 's ring template 
matching at step S1407. 

The calculation section 210 reads the template image from 
the memory 240 in accordance with the size of the specified 

25 Landolt 's ring, and obtains the spatial frequency distribution 
Temp(x, y) thereof (S1301). Next, the calculation section 210 
obtains the two-dimensional Fourier transform FT(u, v) of 
Temp(x, y) (S1303). The calculation section 210 obtains the 
two-dimensional Fourier transform OR(u, v) of the spatial 

30 frequency distribution of the index image data by the simulation 
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of the retinal image, and multiplies OR(u, v) by the spatial 
frequency distribution FT(u, v) of the template as indicated 
by a following expression, and obtains OTmp(u, v) (S1305). 
OR(u, v) X FT(u, v) OTmp(u, v) 
5 The calculation section 210 performs the two-dimensional 

inverse Fourier transform of OTmp(u, v) to obtain TmpIm(X, Y) 
(complex matrix of 4a x 4a) (S1307). The calculation section 
210 acquires the maximum value of the absolute value of TmpIm(X, 
y) to obtain the score n (S1309), 
10 By tcdcing such a correlation, when the simulation index 

image is close to the original image, the score is high, and 
in the case of blur, the score becomes low according to that. 

5-2. Second flowchart (astigmatic correction - 1) of visual 

15 acuity simulation 

Figs. 11 and 12 show second flowcharts (1) and (2) of the 
visual acuity simulation. Figs. 11 and 12 show the flowcharts 
in which the simulation of a retinal image is performed, a 
correction astigmatic eixis A and an astigmatic power C are 

20 obtained so that the Landolt's ring can be detected, and the 
visual acuity at the time of the correction is estimated. This 
example shows a case where the astigmatic power has a negative 
value . 

Similarly to step S1401, the calculation section 210 
25 calculates a tentative spherical power Sr (S1401). By 
specifying the tentative spherical power Sr, in order to avoid 
the retina from approaching the front focal line relatively to 
the rear focal line, it may be set to be slightly weaker 
correction than Sr specified at step S1401 (for example, Sr 
30 Sr + 0.5D) . Alternatively, what is obtained by subtracting the 
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refractive value or 1/2 of the astigmatic power Cs calculated 
from the wavefront aberrations from S obtained as described 
above, or what is specified to be slightly weaker correction 
than them may be used. Next, the calculation section 210 
5 specifies the astigmatic power Cs for simulation (S1501) • For 
example , Cs may be obtained by using the refractive value or 
the astigmatic power C calculated from the wavefront 
aberrations, or a correspondence table storing Cs corresponding 
to correction elements such as S or C or Zernike coefficients 

10 is stored in the memory 240, and it may be obtained by referring 
to that. Next, the calculation section 210 specifies the 
Landolt*s ring of visual acuity Vs (for example, Vs = 0.1) 
similarly to the above (S1453). 

At steps S1405 to S1413, similarly to the above, the 

15 calculation section 210 performs the processing such as the 
Landolt's ring retinal image simulation and Landolt's ring 
template matching. In the case of No at step S1411 or S1413, 
the judgment part 212 judges that the Landolt*s ring can not 
be detected, and stores the direction at this time and that the 

20 detection could not be made in this direction into the memory 
240 (S1419). After step S1419 or in the case of Yes at step 
S1413, similarly to the above, the calculation section 210 
performs the processing of step S1421 and S1455. 

At step S1455, in the case where it is judged that 

25 detection could be made in number of directions equal to or more 
than half of specified number of direction, the calculation 
section 210 stores the specified correction element into the 
memory 240 (S1503). Next, the correction element specifying 
section 213 specifies V = Vs, and specifies the Landolt ' s ring 

30 of visual acuity Vs = Vs + 0.1 (S1505). At this time, with 



respect to the Landolt's ring to be specified, in accordance 
with the specification at step S1453, the Landolt's ring of high 
contrast or low contrast is specified. Thereafter, advance is 
made to S1405, and the image data formation section 211 performs 
5 the retinal image simulation based on the specified correction 
element and the Landolt's ring to obtain index image data, and 
performs the processing subsequent to step S1407. 

On the other hand, in the case of No at step S1455, the 
judgment section 212 judges whether the simulation is performed 

10 in all astigmatic axis angle directions (0 to 180) (S1507). 
Here, in the case of No, the correction element specifying 
section 213 rotates the astigmatic axis angle As (for example. 
As - As + 5) (S1509) • Thereafter, advance is made to step S1453, 
and the processing subsequent to step S1453 is repeatedly 

15 performed. 

Next, referring to Fig. 12, in the case where the judgment 
section 212 makes a judgment of Yes at step S1507, the correction 
element specifying section 213 of the calculation section 210 
substitutes, as the astigmatic axis angle A, As at the time when 

20 the visual acuity V is largest (S1511). Incidentally, when 
there are plural As at the time when it is largest , one at which 
the number of Landolt's rings which could be detected by the 
visual acuity V is largest is specified as the astigmatic eixis 
angle A, and further, when there are plural such As, one is 

25 specified at which the sum of nh in the direction in which 
detection by the visual acuity V could be made becomes maximum. 
By this, the astigmatic axis angle A is determined. 

At steps S1453, and S1405 to S1413, as described in the 
above embodiment, based on the specified Sr, Cs and A, the 

30 calculation section 210 performs the respective processings 
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such as the Landolt ' s ring retinal image simulation and the 
Landolt ' s ring template matching . 

In the case of No at step S1411 or S1413, the judgment 
section 212 judges whether Cs exceeds a previously determined 
5 allowable value (for example, Cs-IOD) (S1515). Here, in the 
case of No, the correction element specifying section 213 
specifies the correction element of Cs to be slightly strong 
(for example, Cs-0.25D) (S1517), the image data formation 
section 211 performs the Landolt 's ring retinal image 

10 simulation based on this correction element (S1405). The 
calculation section 210 repeatedly performs the processing 
subsequent to step S1407 with respect to the index image data 
obtained by this simulation. On the other hand, in the case 
of Yes at step S1515, the judgment section 212 judges that the 

15 Landolt's ring can not be detected (S1419), and stores the 
direction at this time and that the detection could not be made 
in this direction into the memory 240. 

After step S1419 or in the case of Yes at step S1413, 
similarly to the above, the calculation section 210 performs 

20 the processing of step S1421 and S1455. In the case of Yes at 
step S1455, the calculation section 210 performs the processing 
of step S1503 and S1505 . The processing of the respective steps 
is the same as the above. Thereafter, advance is made to step 
S1405, and the image data formation section 211 performs the 

25 retinal image simulation based on the specified correction 
element and the Landolt's ring to obtain the index image data, 
and performs the processing subsequent to step S1407. 

On the other hand, in the case of No at step S1455, the 
calculation section 210 performs the data output (S1423) . That 

30 is, the calculation section 210 displays the visual acuity V 



at this time, the astigmatic power C = Cs, the astigmatic axis 
A, the spherical power S = Sr, the direction in which detection 
could be made, the simulation result and the like on the display 
section 230 and stores them in the memory 240. 

5 

5-3. Third flowchart (astigmatic correction - 2) of visual 
acuity simulation 

Fig. 13 shows a third flowchart of the visual acuity 
simulation. Fig. 13 shows the flowchart in which the MTF is 

10 used as an evaluation parameter to obtain the astigmatic axis 
A and the astigmatic power C, and the visual acuity at the time 
of the correction is estimated. 

In step S1401, as described above, the arithmetic part 
210 calculates a tentative spherical power Sr. Then, the 

15 arithmetic part 210 initially specifies an astigmatic power Cs 
and the angle As of an astigmatic eixis, both of which are 
astigmatic components, and a comparison numeral Mh (S1571). 
These values may be stored in advance in the memory 240, or may 
be input through the input part. The arithmetic part 210 

20 initially specifies, for example, Cs = 0, As « 0, and Mh = 0. 

The arithmetic part 210 calculates the MTF (modulation 
transfer function) (S1573) according to the wavefront 
aberrations obtained before. A specific MTF calculation 
method will be described later. The arithmetic part 210 

25 calculates a comparison niuneral M from an MTF cross -section at 
the specified angle As of the astigmatic cixis (S1575). As the 
comparison numeral M, the total sum of MTF values, an MTF 
cross -section, or the sum of 3, 6, 12, and 18 cpd, for example, 
can be used. The arithmetic part 210 stores currently set As 

30 and M in the memory 240. 



The determination part 212 of the arithmetic part 210 
determines whether M st Mh (S1577). If no, the processing 
proceeds to step S1581. If yes, the correction-factor setting 
part 213 of the arithmetic part 210 sets Mh = M and A = As (S1579) . 
5 Then, the determination part 212 determines whether As Is 180 
or larger (S1581). If no, the correction -factor setting part 
213 rotates the angle As of the astigmatic cixls (for example. 
As = As + 5) (81509). Then, the arithmetic part 210 goes back 
to step S1575 and repeats the processes to obtain the maximum 

10 value of M In an axis -angle range of 0 to 180 degrees and the 
angle As (weak main longitude line or strong main longitude 
line) of the astigmatic axis equal to the direction where the 
maximum value of M Is obtained. 

If yes In step S1581, In other words. If the angle A of 

15 the astigmatic clxIs Is obtained, the arithmetic part 210 

calculates the MTF according to the astigmatic components Cs 
and As = A. The arithmetic part 210 further calculates the 
comparison numeral M from each MTF cross -section at 0 to 180 
degrees (for example, at an Interval of 5 degrees). 

20 The determination part 212 determines whether the 

calculated Ms are almost equal at the angles (S1589). For 
example. It can be determined by determining whether the 
difference between the maximum value of M and the mlnlmiim value 
of M Is smaller than a predetermined threshold. If no In step 

25 S1589, the arithmetic part 210 slightly changes the astigmatic 
power Cs (for example, Cs = Cs - 0.25) (S1591) , and the processes 
of step SI 58 5 and subsequent steps are repeated. If yes In step 
S1589, the arithmetic part 210 specifies C = Cs (S1593). 

The calculation section 210 estimates the visual acuity 

30 after the correction based on the obtained astigmatic power C 



and the astigmatic axis angle A (S1594), For example, the 
calculation section 210 performs the processing of a fourth 
flowchart described later to estimate the visual acuity after 
the correction. Incidentally, instead of estimating the 
5 visual acuity or in addition to the estimation of the visual 
acuity, the calculation section 210 may obtain the contrast 
sensitivity. 

The calculation section 210 stores the obtained 
astigmatic power C and the astigmatic axis angle A into the 
10 memory 240 and displays them on the display section 230 as the 
need arises (S1595), 

5-4. Visual acuity estimation 

Fig. 14 is a fourth flowchart of the visual acuity 

15 simulation of step S107. Besides, the flowchart shown in Fig. 
14 is also a sub-flowchart of the foregoing step S1594. First, 
the calculation section 210 specifies correction data for 
simulation (S1452). For example, the calculation section 210 
can use, as the correction data, a value calculated based on 

20 the refractive value or wavef ront aberrations , or the spherical 
power S obtained as stated above, the astigmatic power C, or 
the astigmatic axis angle A. Besides, by setting the respective 
elements of the correction data to be 0, the calculation section 
210 can estimate the visual acuity in the environment of the 

25 subjective eye when the correction is not made. Besides, for 
example, the calculation section 210 may specify the astigmatic 
power C, the astigmatic axis angle A, and/or the spherical power 
correction S, which are obtained by performing the processing 
of the foregoing flowchart. Since the processing of each step 

30 subsequent to step S1453 is the same as the processing denoted 
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by the same reference character in the flowcharts shown in Figs . 
7 and 11, the detailed description will be omitted. 

5-5. Contrast sensitivity 
5 The calculation section 210 can calculate contrast 

sensitivity as the visual acuity simulation of step S107. The 
calculation section 210 obtains Mopt(r, s) of the eye optical 
system based on the wavefront aberrations, and calculates the 
contrast sensitivity based on the obtained MTF. Besides, the 

10 calculation section 210 displays the calculated contrast 
sensitivity on the display section 230 or stores it into the 
memory 240. Incidentally, the contrast sensitivity is not only 
calculated in the processing of step S107, but also can be 
calculated during the processing of the first to fourth 

15 flowcharts and can be displayed. 
(MTF calculation) 

First, the calculation of the MTF (Modulation transfer 
function) will be described. 

The MTF is an index indicating the transfer 

20 characteristic of the spatial frequency, and is widely used to 
express the performance of an optical system. In this MTF, 
visibility can be predicted by obtaining, for example, the 
transfer characteristic for 0 to 100 sinusoidal gray lattices 
per degree. In this embodiment, as described below, a 

25 monochromatic MTF may be used or a white MTF may be used. 

First, the monochromatic MTF is calculated from the 
wavefront aberration W(x, y) . Incidentally, W(x, y) is an input 
value (measured value), and with respect to the corneal 
aberrations, the corneal wavefront aberrations obtained from 

30 the cornea shape can also be used. 
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When the monochromatic MTF is obtained, the calculation 
section 210 obtains the pupil function f(x, y) from the 
wavefront aberrations as described below. 



5 (i: imaginaory number, k: wave number vector (2jt/X), X: 
wavelength) 

At this time, in view of the Stiles -Crawford effect, (e"^)^ 
(ais, for example, approximately 0 . 05) may be multiplied. Here, 
rp is the radius of the pupil. 
10 The calculation section 210 performs the Fourier 

transform of this pupil function f (x, y) to obtain an amplitude 
distribution U(u, v) of point images as indicated by a following 
expression. 



15 ( X : wavelength 

R: distance from the pupil to an image point (retina) 

(u, v) : coordinate value in a plane perpendicular to the optical 

axis while the image point O is made the origin 

(x, y) : coordinate value in the pupil plane) 

20 The calculation section 210 multiplies U(u, v) by its complex 
conjugate and obtains I(u, v) as the point spread function (PSF) 
by a following expression. 

I(u, V) = U(u, v)U*(u, V) 
Next, the calculation section 210 performs the Fourier 

25 transform ( or autocorrelation) of PSF to perform normalization, 
and obtains OTF. 



f (X, y) = e^'^*^^' y> 



00 



j_2£ 
R Jl 



Amplitude U(u,v)= JJf(x,y)exp 



(ux + vy) dxdy (14) 



-00 
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R(r,s) = J Jl(u,v)e-'2.^™"^Mudv 

-00 

(r, s: variables in spatial frequency area) (15) 

~ R(0,0) 

Besides, since the magnitude of the OTF is the MTF, 
MTF(r, s) = |OTF(u, v)| 
is established. 

5 The white-color MTF is calculated from the single-color 

MTF, obtained as described above. 

To obtain the white-color MTF, the MTF is weighted at each 
wavelength and added. Since the above -described MTF has a 
different value at each wavelength, the MTF can be expressed 
10 in the following way when the MTF at a wavelength X is indicated 
by MTFx. 



MTF(r,s) = 



(16) 

The MTF is highly weighted at visible- light wavelengths, 
15 and the calculation is made. 

More specifically, the MTF is obtained in the following 
way when it is assumed, for example, that the three primary 
colors (R, G, and B) are specified such that red light has a 
wavelength of 656.27 nm with a weight of 1, green light has a 
20 wavelength of 587.56 nm with a weight of 2, and blue light has 
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10 



15 



a wavelength of 486.13 nm with a weight of 1. 

MTF(r, s) = (1 X MTF656.27 + 2 x MTF587.56 + 1 x MTF486.13) / ( 1 
+ 2 + 1) 

Since the white -light MTF is measured only at one 
wavelength (840 nm) , calibration may be performed for other 
wavelengths according to the result of measurement, as 
compensation, to obtain the MTF at each wavelength. More 
specifically, when the eye optical characteristic measuring 
apparatus measures eye aberration, for example, at 840 nm, color 
aberration Wa(x, y) corresponding to a shift from the wavefront 
aberrations W84o(x, y) at a wavelength of 840 nm is measured with 
the use of an eye model, W84o(x, y) is added to the color 
aberration Wa(x, y ) , and the MTF is calculated at each wavelength 
from this wavefront aberrations in the following way. 

Wx(x, y) := W84o(x, y) + Wa(x, y) 



20 



(Contrast sensitivity calculation) 

Next, the contrast sensitivity will be described. The 
contrast sensitivity is expressed by a following expression 
(Peter G.J. Barten, "Contrast Sensitivity of the Human Eye and 
Effects on Image Quality" , SPIE Optical Engineering Press 
1999) . 



S{r,s) = 




(17) 



25 Where, the respective pareuneters are as follows: 
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Mopt{r# s): MTF of the eye optical system, k: S/N ratio:3, T: 
weighting time of a nervous system: 0.1 sec, X©: visual angle 
of a matter: 3.8 deg, Xmax: maximum visual angle of space 
weighting: 12 deg, Nmax: maximum frequency when weighting Is 
5 performed: 15 cycles, Ti: quantum efficiency of a light receptor 
of an eye: 0.3, p: photon conversion coefficient (ORT) of a light 
source: 1.24 (liquid crystal Is also acceptable), E: retina 
Illumination (Troland) : 50 (cdm^) x r^Ji (mm) = 50r^jt(td) (r: pupil 
radius) 100 or less, <|)o: spectral density of nervous system 

10 noise: 3 x 108 sec*deg^, Uq: spatial frequency of side 
suppression: 7 cycles/deg. By using this expression, not the 
contrast sensitivity by the eye optical system, but the contrast 
sensitivity of the whole optic system including other elements 
(for example, nervous system) can be predicted. 

15 Fig. 15 is an explanatory view of the contrast sensitivity. 

In the graph shown in Fig. 15, the vertical axis indicates the 
contrast sensitivity calculated by using the foregoing 
expression, and the horizontal axis indicates the spatial 
frequency, and the graph is a one- dimensional graph (graph at 

20 the time of, for example, s = 0) on a certain section passing 
through the origin. By obtaining the contrast sensitivity of 
the whole optic system corresponding to the spatial frequency, 
for example, the visibility of a stripe index can be predicted. 
Besides, an eye doctor or the like can compare, for 

25 example, the contrast sensitivity displayed on the display 
section with the sensitivity by subjective measurement. For 
example, it is possible to compare the sensitivities of 3cpd, 
6cpr, 9cpd and 12cpd in the x direction obtained by general 
subjective measurement and according to vertical stripe indexes 

30 with the contrast sensitivities corresponding to the respective 
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spatial frequencies at the time of s = 0, Incidentally, since 
the contrast sensitivity does not depend on the angle in the 
case where it is displayed in polar coordinates and is rotation 
symmetry, it can also be displayed while the horizontal axis 
5 of the graph is made the amplitude component of the polar 
coordinate display. 

Incidentally, the foregoing first, second, third, and 
fourth flowcharts and the calculation of the contrast 
sensitivity are combined and are used to obtain the correction 

10 values of the spherical power, the astigmatic power, and the 
astigmatic axis, and the visual acuity and/or the sensitivity 
at the time of the correction may be measured. In the case where 
the astigmatism is considered, since the spherical power 
calculated by the first flowchart becomes an equivalent 

15 spherical power Se, the spherical power is made S = Se - (1/2)C. 

6. Correction image simulation 

6-1 First flowchart of correction image simulation 

Fig. 23 is a flowchart of the correction image simulation 

20 of the step S107. 

The arithmetic part 210 calculates a best image condition 
(S201). As described later, the details are such that the 
arithmetic part 210 obtains a lower order Zernike coefficient 
so that the Strehl ratio becomes msuclmum or the phase shift 

25 becomes as small as possible, and obtains corrective correction 
data. As the corrective correction data, suitable data can be 
named eimong, for example, coefficients corresponding to defocus, 
astigmatism components, S, C, A, higher order spherical 
aberrations, higher order astigmatism aberrations, higher 

30 order coma aberrations, the Strehl ratio and the like. 
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The arithmetic part 210 obtains the wavef ront aberrations 
W(x, y) at the time of the best image condition, and calculates 
the pupil function f(x, y) from W{x, y) by the following 
expression (S203). Detail description of processes from step 
5 S205 to S211 is omitted, because each steps are same or similar 
to same number of steps of Fig. 8. 

The arithmetic part 210 displays the LandImage(X, Y) and 
PSF(X, Y) on the display part 230 by a suitable display method 
of a drawing, graphic data, a graph and/or a numerical value, 
10 and suitably stores the data in the memory 240 (S213). The 
arithmetic part 210 reads out corrective correction data from 
the memory 240 as the need arises, and outputs it to the display 
part 230 (S215). 

15 6-2 Calculation of correction data based on Strehl ratio 

Fig. 24 shows a flowchart concerning a first example of 
best image condition calculation. Fig. 24 is the detailed 
flowchart concerning the foregoing step S201. 

First, the arithmetic part 210 sets a threshold value for 

20 respective aberration quantities RMSi^^"^ as a branch condition 
(S401). For example, this threshold value can be made a 
sufficiently small value (for example, 0.1) of aberration. The 
arithmetic part 210 calculates the Zernike coefficients Ci^^"^ 
from the measured detection wavef ront, and converts them to the 

25 aberration quantities RMSi^^"^ by the following expression 
(S403) . 




i 

2(/ + 1) 



2j-i 



(ef'-' = 2 (2 J = i) , sf'-' = 1 (2 J * i) ) 



(18) 



The arithmetic part 210 judges whether at least one of 
the values of RMSi^^"^ (i > 2) is the threshold value or higher 
(S405). Here, in the case where a judgment of No is made, it 
proceeds to step S419. On the other hand, here, when a judgment 
of Yes is made, the arithmetic part 210 carries out a next 
processing. 

That is, the arithmetic part 210 judges whether at least 
one of the higher order spherical aberration quantities R4^, 
Re® — of aberration quantities RMS (Ri^^"^) is the threshold value 
or higher (S407) . Here, in the case of Yes, the arithmetic part 
210 causes the aberration to change a coefficient (C2®) 
corresponding to the def ocus so that the Strehl ratio becomes 
maximum (S409), and on the other hand, in the case of No, it 
proceeds to step S411. Next, the arithmetic part 210 judges 
whether at least one of the asymmetrical higher order coma 
aberration quantities RMSi^^"^ (i: odd number) is the threshold 
value or higher (S411) • Here, in the case of Yes, the arithmetic 
part 210 causes the aberration to change the coefficient (C2®) 
corresponding to the def ocus so that the Strehl ratio becomes 
maximum (S413), and on the other hand, in the case of No, it 
proceeds to step S415. Further, the arithmetic part 210 judges 
whether at least one of the higher order astigmatism aberration 
quantities RMSi^^"^ (i: even number and 2j-l 0) is the threshold 
value or higher (S415) • Here, in the case of Yes, the arithmetic 
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part 210 adds astigmatism components (Ca"^, Cz^) to the 
aberration so that the Strehl ratio becomes maximum (S417) , and 
on the other hand, in the case of No, it proceeds to step S419, 
In this way, the arithmetic part 210 calculates OTF(u, 
5 v) and PSF(X, Y) from the aberrations, and further calculates 
the corrective correction data (suitable data such as 
coefficients corresponding to the defocus, astigmatism 
components, S, C, A, higher order spherical aberrations, higher 
order astigmatism aberrations, higher order coma aberrations, 

10 and Strehl ratio) from the Zernike coefficients, and stores them 
in the memory 240 (S419), 

Incidentally, in order to correct only a desired 
component among the defocus and the astigmatism components, any 
of the pairs of the steps S407 and S409, the steps S411 and S413, 

15 and the steps of S415 and S417 may be omitted, or a step may 
be added to correct suitable higher order aberrations or Zernike 
coefficients other than these. For example, in the case where 
a fourth- order spherical aberration is mainly included in the 
higher order aberrations, the corrective correction data can 

20 be obtained by correcting in the direction in which the defocus 
quantity corresponding to the lower order aberrations are 
increased. 

Next, the detailed processing of the steps S409, S413 and 
S417 will be described. In the respective steps, the arithmetic 
25 part 210 carries out the processing as follows. 

In order to obtain a more suitable image plane, the 
arithmetic part 210 adds to the wavefront aberrations W(x, y) 
the lower order Zernike coefficients Ci^^"^(l ^ ± ^ 2) at each 
step presently noted for the aberration quantities comparable 
30 to the higher order aberration quantities according to the 



threshold value of the last noted higher order aberration 
quantities (RMS4°, RMSe^ - ) in the flow. For example, €2° is 
added at the step S409; Cz^ , at the step 413; and €2"^^ Ca^. at 
the step S417. 

Further, the pupil function f (x, y) is obtained from the 
wavefront aberrations in the manner described below. 

f (x, y) = e^^*^'^' y> 

(i: imaginary number, k: wave number vector (2jt/X), X: 
wavelength) 

The arithmetic part 210 performs the Fourier 
transformation on this pupil function f(x, y) , so that an 
amplitude distribution U(u, v) of a point image is obtained as 
in the following expression. 

amplitude u(u,v) = j]f(x,y)exp|^-l^(ux H- vy)j dxdy 

(19) 

(X: wavelength 

R: a distance from a pupil to an image point (retina) 
(u, v) : a coordinate value on a plane orthogonal to an optical 
axis while an image point O is made the origin 
(x, y) : a coordinate value on a pupil plane) 

The arithmetic part 210 multiplies U(u, v) by its complex 
conjugate and obtains I(u, v) as a point image intensity 
distribution (PSF) by the following expression. 

I(u, v) = U{u, v)U*(u, V) 

Besides , when the center intensity of PSF at the time when 
there is no aberrations (W(x, y) = 0) is made Io(0, 0) , the Strehl 
ratio is defined as follows: 

Strehl ratio = 1(0, 0)/Io(0, 0). 
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In the first example, the arithmetic part 210 recursively 
or analytically obtains a value of the lower order Zernike 
coefficient C13 (l^i^2) so that the value of the Strehl ratio 
becomes meLximum. 

5 

6-3 Calculation of correction data based on a phase shift 

Next, Fig. 25 is a flowchart concerning the second example 
of the best image condition calculation. 

First, the arithmetic part 210 sets a threshold value for 

10 the respective aberration quantities RMSi^^'*^ as a branch 
condition (S501). For example, this threshold value is made a 
sufficiently small value (for example, 0.1) of aberration. 

The arithmetic part 210 calculates the Zernike 
coefficients Ci^^'^ from the measured detection wavefront, and 

15 converts them to the aberration quantities RMSi^^**^ by the 
expression indicated in the first example (S503). The 
arithmetic part 210 judges whether at least one of RMS i^^'^ (I 
> 2) is the threshold value or higher (S505) . Here, in the case 
of the judgment of No, it proceeds to step S519. On the other 

20 hand, in the case of the judgment of Yes, the arithmetic part 
210 carries out a next processing. 

That is, the arithmetic part 210 judges whether at least 
one of the higher order spherical aberration quantities R4°, 
Re^ — is the threshold value or higher (S507) . Here, in the case 

25 of Yes, the arithmetic part 210 causes the aberration to change 
the coefficient (€2°) corresponding to the defocus so that the 
phase shift becomes as small as possible (S509 ) , and on the other 
hand, in the case of No, it proceeds to step S511. Next, the 
arithmetic part 210 judges whether at least one of the higher 

30 order coma aberration quantities RMSi^^"^ (i: odd number) is the 
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threshold value or higher (S511) . Here, in the case of Yes, the 
arithmetic part 210 causes the aberration to change the 
coefficient (Ca^) corresponding to the defocus so that the phase 
shift becomes as small as possible (S513) , and on the other hand, 
5 in the case of No, it proceeds to step S515. Further, the 
arithmetic part 210 judges whether at least one of the higher 
order astigmatism aberration quantities RMSi^^"^ (i: even number 

and 3 5^ 0) is the threshold value or higher (S515). Here, in 
the case of Yes, the arithmetic part 210 adds the astigmatism 

10 components (Ca"^* Cz^) to the aberration so that the Strehl ratio 
becomes meiximum (S517), and on the other hand, in the case of 
No, it proceeds to step S519. 

In this way, the arithmetic part 210 calculates OTF(u, 
v) and PSF(X, Y) from the aberrations, and further calculates 

15 the corrective correction data (suitable data such as the 
coefficients corresponding to the defocus, astigmatism 
components, S, C, A, higher order spherical aberrations, higher 
order astigmatism aberrations, higher order coma aberrations, 
and Strehl ratio) from the Zernike coefficients , and stores them 

20 in the memory 240 (S519). 

Incidentally, any of the pairs of the steps S507 and S509, 
the steps S511 and S513, and the steps S515 and S517 may be 
omitted so that only a desired component is corrected among the 
defocus and the astigmatism components. Besides, a step may be 

25 added so that suitable higher order aberrations or Zernike 
coefficients other than these is corrected. 

Next, the detailed processing of the steps S509, S513 and 
S517 will be described. The arithmetic part 210 carries out the 
processing as follows. 

30 First, as described in the detailed processing of the 
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steps S409, S413 and S417, the arithmetic part 210 obtains the 
point image intensity distribution (PSF) from the expression 
of the wavefront at the time of the objective complete 
correction calculated from the Zernlke coefficients • Next , the 
arithmetic part 210 performs a Fourier transformation (or 
autocorrelation) on the PSF to normalize it as in the following 
expression and obtains OTF. 

oo 

R(r,s) = J Jl(u,v)e-'2'^~*'"Mudv 

— oo 

(r, s: a variable of a spatial frequency region) 

R(r,s) 



OTF = 



R(0,0) 



(20) 

10 In general, the amplitude of a spatial frequency region 

and a phase distribution R(r, s) become complex numbers, and 
when its real number part is A(r, s), and its imaginary part 
is B(r, s) , 

R(r, s) = A(r, s) + iB(r, s) 

15 and the shift of the phase (phase shift, PTF) is 

calculated by 

<?>(r,j) = tan — 

(21) 

In the second example, the arithmetic part 210 
20 recursively and analytically obtains such a value of the lower 
order Zernike coefficient Ci^^"^ that a value at which the R(r, 
s) has an extreme value is brought to a high frequency to the 
extent possible, that is, the phase shift becomes as small as 
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possible. 

Incidentally, with respect to the first example and the 
second example of the best Image condition calculation, both 
the processings may be carried out to obtain such a condition 
5 that the Strehl ratio Is large, and the phase shift Is small. 

7. Display Excunple 

?•! Display example at the time of visual acuity measurement 
Fig. 16 shows, with respect to the best Image display - 

10 template matching and as numerical data, the spherical power 
S, astigmatic power C, and astigmatic axis angle Ax, which are 
compensation correction data, the pupil dlcuneter, and the 
spherical power S, astigmatic power C, astigmatic axis euigle 
Ax, and corrected visual acuity, which are measured values 

15 before the compensation correction. Further, corrected visual 
acuity In the correction data may be displayed. In this example, 
since the component of the higher order aberration has a 
predetermined value or more, a difference In numerical value 
occurs between the compensation correction data and the 

20 measured values • In these drawings, the wavefront aberration, 
PSF, OTF, OTF (two-dimensional display), S, C, Ax, Landolt's 
ring, visibility of the Index and the like are displayed on the 
display section 230. Further, for example, the contrast 
sensitivity shown In Fig. 15 may be displayed on the display 

25 section 230. Besides, some of them may be suitably selected 
and displayed. 

Fig. 17 Is a view showing a display example concerning 
the comparison between pre-compensatlon and post -compensation. 
This drawing shows the corrected visual acuity, wavefront 

30 aberration, and visibility of the Landolt's ring, which are 



obtained before and after the compensation, and the pupil 
diameter. As shown in the drawing, it is indicated that after 
the compensation correction, the wavefront aberration becomes 
relatively uniform, and the Landolt's ring is relatively well 
5 seen. Besides, the corrected visual acuity of the subjective 
eye under an environment after the compensation correction is 
shown. 

Fig. 18 is an explanatory view of an example of 
prescription data for eyeglasses/contacts. Fig. 19 is an 

10 explanatory view of an example of data for refractive surgery. 

The respective data are stored in the memory 240 from the 
calculation section 210, and/or displayed on the display 
section 230. This example indicates that in the data of the 
case where the refractive surgery is performed while only the 

15 SCA is made the compensation correction data, the corrected 
visual acuity is improved by performing the correction in such 
a way that the value of S in the compensation correction data 
is intensified, the value of C is weakened, and the axial 
direction of A is slightly changed. Besides, Fig. 19 shows 

20 expected values of the respective parameters under the 
illustrated pupil diameter after the compensation correction. 

Fig. 20 is an explanatory view of an example of 
prescription data for eyeglasses /contacts when the 
environmental condition is changed. For example, the pupil 

25 diameter of the eye 60 to be measured is measured in the 
illumination states corresponding to the respective 
environmental conditions , and the correction data and corrected 
visual acuities at the respective pupil diameters are displayed. 
It is indicated that the compensation correction data slightly 

30 varies according to the pupil diameter. That is, it is 



indicated that the optimum prescription value varies according 
to the environment of the subjective eye. Besides , for example, 
the doctor considers the environment of the subjective eye, and 
can select the prescription value. Incidentally, the 
5 environmental conditions to be displayed can be suitably 
changed • 

In the example shown in Fig. 20, although the correction 
data corresponding to the respective environmental conditions 
are obtained, and the visual acuities under the environmental 

10 conditions are displayed, the visual acuity under another 
environmental condition can also be estimated and displayed. 
For example, in the case where correction is made by the 
compensation correction data in the daytime, the visual 
acuities in the daytime, under a fluorescent lamp and in a room 

15 and daytime can also be predicted and displayed. 

Fig. 21 is an explanatory view of an example of pupil data 
when the environmental condition is changed. For example, the 
pupil diameter of the eye 60 to be measured is measured under 
illumination states corresponding to the respective 

20 environmental conditions, and the shift amount (x direction, 
y direction) from the limbus center of the pupil center at the 
respective pupil diameters and the corrected visual acuity are 
shown. It is indicated that the pupil center is shifted by the 
change of the environmental condition, and the center (origin) 

25 at the time of analysis is shifted. 

Fig. 22 is a comparison view of prescription data for 
eyeglasses /contacts and measurement at constant pupil 
diameters. There are shown, for example, the correction data 
and corrected visual acuities in the case where similarly to 

30 the conventional measurement, the pupil diameter is made 4mm 



and 6mm and the case where the pupil diameter is measured (for 
example, under illumination of 50 Ix) . The correction data and 
corrected visual acuity slightly vary between the case where 
the pupil diameter is fixed and the case where it is measured. 
5 Incidentally, although Fig. 21 shows, as an example, data using 
the pupil diameter under the illuminated of 50 Ix, it is possible 
to estimate the visual acuity of the subjective eye under an 
appropriate environment by suitably changing the illumination 
condition. 

10 Incidentally, in the foregoing drawing, although the 

visual acuity is expressed by the decimal visual acuity, it may 
be displayed by the logMAR visual acuity. Besides, the 
condition to be displayed can be suitably changed. 

15 7.2 Display example in correction data simulation 

Fig, 26 shows, with respect to the best image display - 
Strehl optimization and as numerical data, the spherical power 
S, astigmatic power C, and astigmatic axis angle Ax, which are 
compensation correction data, the spherical power S, astigmatic 

20 power C, and astigmatic 3lx±b angle Ax, which are measured values 
before the compensation correction, and the pupil dicuneter used 
for the calculation of the correction data. In this example, 
since the components of the higher order aberrations have 
predetermined values or more, a difference in numerical value 

25 occurs between the compensation correction data and the 
measured values. 

Fig. 27 shows, with respect to the best image display - 
PTF optimization and as numerical data, the spherical power S, 
astigmatic power C, and astigmatic axis angle Ax, which are 

30 compensation correction data, the spherical power S, astigmatic 
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power C, and astigmatic axis angle Ax, which are measured values 
before the compensation correction, and the pupil diameter used 
for the calculation of the correction data. In this example, 
since the components of the higher order aberrations have 
5 predetermined values or more, a difference in numerical value 
occurs between the compensation correction data and the 
measured values . 

Fig. 28 Is a view showing a display example concerning 
the comparison between pre -compensation and post -compensation. 

10 This drawing shows the wavef ront aberrations before and after 
the compensation, the visibilities of the Landolt's ring, 
Strehl ratios, and the pupil diameter. As shown in the drawing, 
the Strehl ratio after the compensation is higher, the wavef ront 
aberrations become relatively uniform, and the Landolt's ring 

15 can also be relatively well seen. 

Fig. 29 is an explanatory view of an example of 
prescription data for eyeglasses/contacts. Fig. 30 is an 
explanatory view of an example of data for refractive surgery. 
The respective data are stored from the calculation 

20 section 210 into the memory 240, and/ or are displayed on the 
display section 230. This exeunple indicates that in the data 
of the case where the refractive surgery is performed while only 
the SCA are made the compensation correction data, the Strehl 
ratio becomes high and the correction effect becomes high by 

25 performing the correction in such a manner that the value of 
S in the compensation correction data is intensified, the value 
of C is weeikened, and the axial direction of A is slightly 
changed • 

Fig. 31 is an explanatory view of an excunple of 
30 prescription data for eyeglasses/contacts when the 
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environmental condition is changed. For example, the pupil 
diameters of the eye 60 to be measured are measured under 
illumination states corresponding to the respective 
environmental conditions, and the correction data at the 
5 respective pupil dieuneters are displayed. It is indicated that 
the compensation correction data slightly varies according to 
the pupil diameters. That is, it is indicated that the optimum 
prescription value varies according to the environment of the 
subjective eye. Besides, for example, the doctor or the like 

10 considers the environment of the subjective eye, and can select 
the prescription value. Incidentally, the environmental 
condition to be displayed can be suitably changed. 

Fig. 32 is a comparison view of prescription data for 
eyeglasses /contacts and measurement at constant pupil 

15 diameters. For example, correction data when the pupil 
diameter is made 4 mm and 6 mm similarly to the conventional 
measurement, and correction data of the case where the pupil 
diameter is measured (for example, under illumination of 50 Ix) 
are displayed. The correction data slightly varies between the 

20 case where the pupil diameter is fixed and the case where it 
is measured. In this embodiment, it is possible to obtain the 
optimum correction data corresponding to the pupil diameter of 
the subjective eye. Incidentally, the condition to be 
displayed can be suitably changed. 

25 

8. Modified example 

A modified example of the invention will be described 

below. 

This modified exeunple modifies the calculation method of 
30 the best image condition at S201 of Fig. 23. 
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A component of an 1-th row and a J-th column of Jacoblan 
matrix A Is 

A = > 

Where, fi(x) Is the Strehl ratio, the PTF corresponding 
5 to a suitable frequency, or some values of the PTF corresponding 
to plural frequencies. Besides, It may be a combination of the 
Strehl ratio and the PTF. Besides, a vector x Is an adjustable 
parameter, and here, the sphere (or defocus corresponding to 
that) and two astigmatisms correspond to that. 
10 The calculation expressions of the Strehl ratio and the 

PTF are already given. The Ideal value of the Strehl ratio Is 
1. It Is assumed that the following expression expresses the 
Strehl ratio. 

fl(x) = fs{hC, C2^ C2■^ C2^) (23) 

15 Where, fx denotes the expression of the same Indication 

In the expression (22). 

Besides, for excunple, as the PTF, values corresponding 
to the spatial frequency of 3 cpd, 6cpd, 12 cpd, and 18 cpd are 
taken, and It Is Ideal that this Is 0. 

20 f2(x) = fpTF3(hC, C2°, C2'^ Cz^ ) (24) 

f3(x) = fpTF6(hC, C2°, Cz'^ Cz^ ) (25) 

f4(x) = fpTF12(hC, C2^, C2■^ Cj^ ) (26) 

f5(x) = fpTF18{hC, Cz"", C2"^, C2^ ) (27) 

In the expressions (23), (24), (25), (26) and (27), he 
25 denotes a vector of higher order aberration coefficients, 

denotes a coefficient of a defocus term relating to the sphere, 
C2"^ and C2^ denote coefficients of terms relating to astigmatism. 
The vector he Is given by wavefront aberration measurement, and 
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here. It is constant. Thus, the remaining three coefficients 
are made a parameter vector x and are suitably moved to guide 
fg to the minimum value, which is a task here. 

Here, the partial differentiation of the expression (22) 
5 can be calculated by slightly moving the parameters to prepare 
a change table, and the Jacobian matrix in this system is 
obtained. 

Now, when the task here is expressed in other words again, 
since nonlinear optimization in the case where the Jacobian, 
10 that is, the partial differential coefficient is known has only 
to be performed, when optimizing algorism of a Newton method 
system is used, it is easy to obtain a solution since the example 
is simple. When a specific solution according to a corrected 

Marguardt method is stated, a correction vector Ax can be 
15 obtained by 

(A^WA + XI) Ax = A*W(y - f(x)) (28) 
Here, t at the shoulder of the matrix denotes a transposed 
matrix, and W denotes a weighting matrix. The first element of 
y corresponds to Strehl ratio, and the remainder corresponds 
20 to four components of PTF, it has only to be made (1,0,0,0,0)*. 
X is called a damping factor, and it is made large at first, 
and then, it is made small in accordance with going of 
optimization . 
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(29) 



25 A subscript corresponds to a subscript of f . Weighting 
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suitable for the object of a prescription can be freely 
performed, for example, when the Strehl ratio Is desired to be 
selectively optimized, Wi Is made large. The expression (7) Is 
applied several times, and when 
5 S = W(y - f(x)) (30) 

becomes suitably small (when a conversion condition Is 
satisfied) , the calculation Is stopped, and x at that time Is 
made the solution. By this, the optimum the sphere (or defocus 
corresponding to that) and two astigmatisms are obtained. 

10 Incidentally, the best Image condition calculation can 

also be performed by finding out the position where the Strehl 
ratio becomes maximum or the phase shift (PTF) becomes 
substantially zero while the defocus amount and/or the 
astigmatism component is changed slightly. Besides, a position 

15 where the Strehl ratio becomes meucimum or the phase shift (PTF) 
becomes substantially zero may be obtained by using a well-known 
Newton method. 

9 . Appendix 

20 The apparatus and system of the ophthalmic data 

measurement according to the invention can be provided through 
an ophthalmic data measurement program for causing a computer 
to execute the respective procedures, a computer readable 
recording medium on which the ophthalmic data measurement 

25 program is recorded, a program product including the ophthalmic 
data measurement program and loadable into the internal memory 
of a computer, a computer, such as a server. Including the 
program, or the like. 

Besides, although the measurement data indicating the 

30 refractive power distribution of the eye to be measured is 
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obtained by the optical system 100 shown in Fig, 1 , no limitation 
is made to this, and the structure can be made by another 
aberrometer or the like. 



5 Industrial Applicability 

According to the invention, the optical characteristic 
corresponding to the pupil diameter of the eye to be examined 
and the correction data close to the optimum prescription value 
are calculated, and more accurate measurement can be performed. 

10 Besides, according to the invention, in the results of 

measurement by the eye characteristic measuring appeiratus which 
can measure higher order aberrations, in the case where the 
higher order aberrations are included, the lower order 
aberrations corresponding to the time of the objective complete 

15 correction are not made the compensation correction data, and 
for example, the optical performance is evaluated by the Strehl 
ratio or the phase shift, the lower order aberration amount to 
increase the Strehl ratio and/or to reduce the phase shift is 
calculated, and the compensation correction data, such as S, 

20 C, A, at that time are obtained, so that the correction data 
close to the optimum prescription value of eyeglasses/contacts 
can be obtained. Further, the simulation of the visibility of 
the index for eye examination is performed, and the appropriate 
correction element is obtained, so that it is also possible to 

25 obtain the correction data close to the subjective value. 

According to the invention, the visual acuity of the eye 
to be exeunined can be estimated under the brightness (for 
example, in the daytime or in a room) corresponding to the 
environment of the subjective eye in daily life. Besides, 

30 according to the invention, in view of the pupil diameter of 
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the eye to be examined in daily life, the visual acuity with 
respect to the index of high contrast and/or low contrast can 
be estimated. According to the invention, the contrast 
sensitivity in view of the pupil diameter can be predicted. 
Further, according to the invention, by using the pupil diameter 
under the brightness (for example, in the daytime or in a room) 
corresponding to the environment of the subjective eye, the 
correction data close to the optimum prescription value under 
the environment is obtained, and the visual acuity under the 
environment of the subjective eye can be estimated by the 
correction using the obtained correction data. Besides, the 
simulation of an index such as the Landolt ' s ring on the retina 
in view of the size of the pupil area calculated in the halfway 
process is also singly effective. 



69 



